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Abstract 

Microcanonical equations for several thermodynamic properties of a system, suitable for molec- 
ular dynamics simulations, are derived from the nonextensive Tsallis entropy functional. Two 
possible definitions of temperature, the usual one and a "physical" modification which satisfies the 
zeroth law of thermodynamics, are considered, and the results from both choices are compared. 
Results for the ideal gas using the first definition of temperature are provided and discussed in re- 
lation with the canonical results reported in the literature. The second choice leaves most formulae 
unchanged from their extensive (Shannon-Boltzmann-Gibbs) form. 



PACS numbers: 05.10.-a 



I. INTRODUCTION 



In the last two decades there has been a great deal of interest in nonextensive entropies 
to explain physical phenomena such as anomalous diffusion, believed to be outside the scope 
of the conventional and highly successful Shannon-Boltzmann-Gibbs (SBG) entropy [1]. 
Among the various reported nonextensive functionals, the Tsallis entropy |2j has received 
considerable attention. It is defined as 



S q{Pi} = -kB^PihgqPi, (la) 
i 

where 

log x = lim — — , (lb) 

q'->q 1 — q 

Pi being the probability of each microstate I accessible to the system and ks the Boltzmann 
constant. This definition can be considered an uniparametric generalization of the SBG 
functional, which is recovered when the entropic index q equals one, as can be easily verified 
using L'Hopital rule. A further generalization, not studied in this paper, replaces ks with a 
generic k (q) under the weaker constraint k (1) = ks- 

Most discussions about the Tsallis entropy and the theoretical and computational results 
from its aplication have been made in the context of the canonical ensemble, although a sat- 
isfactory definition of the mean energy has taken some effort to achieve {2]. As with ordinary 
extensive statistics, this ensemble is often more amenable to theoretical calculations. How- 
ever, the microcanonical ensemble is clearly more directly accessible to molecular dynamics 
(MD) simulation, which consists simply in the integration of Newton's Second Law, and as 
such gives rise automatically to the conservation of energy, but lacks the intrinsic notion of 
a thermostat. Thus, it is necessary to have a formalism which allows for the obtention of 
the thermodynamic properties of a system from microcanonical averages. 

Several methods have been devised in order to obtain these formulae for the SBG func- 
tional, adjusted to different needs. For instance, Ray and Graben [3|, guided by didactic 
concerns, use a method based on fluctuations which involves some arbitrary definitions and 
is only valid for relatively large systems. A more systematic and general method, based on 
Laplace transforms, was put forward by Pearson, Halocioglu and Tiller [4j. However, the use 



2 



of this transform is only a method of integration. The present paper applies an equivalent 
but more straightforward method to the Tsallis entropy, and it is structured as follows: in 
the next section, after reviewing some well-know features of the microcanonical ensemble, a 
set of formulae for thermodynamic properties such as the heat capacities and compressibility 
coefficient is developed using the conventional definition of temperature. These formulae are 
applied to the nonextensive ideal gas in order to obtain its thermal and caloric equations. 
After discussing an important shortcoming of this definition of temperature, the previously 
developed results are changed in order to reflect a physical definition of temperature. Finally, 
the main conclusions are summarized. 



II. FORMULAE OBTAINED USING THE ORDINARY TEMPERATURE 

The system to be studied consists of iV classical particles of mass m, with coordinates 
(f N ,p N ) in 6iV-dimensional phase space, in a recipient of fixed volume V and with a to- 
tal energy E, resulting from a Hamiltonian which must be separable into a kinetic and a 
potential part which depend only on the positions and momenta, respectively: 

E = H (f^) = E c (p N ) + U (f N ) . (2) 

The probability distribution over the microstates compatible with these restrictions can be 
obtained using the maximum entropy (MaxEnt) method, by maximization of ( JTatl submitted 
to the normalization ^2 t Pi = 1. In this way, it is trivially obtained that all the microstates 
compatible with the specified thermodynamic coordinates are equiprobable. In practice, 
the energy of the system can only be determined (and constrained) to be in an interval 
[E — SE, E], with SE <C E for a useful measurement. Thus, the number of microstates 
available to the system would be T (E, V, N) — T (E — SE, V, N), with: 

T{E,V,N) = ^ je[E-H (r N ,p N )] dr N df. (3) 

The spatial limits of integration are V N . The momenta can be integrated over the whole 
M. 3N . 9 is the Heaviside step function which restricts the integration to the volume in phase 
space where the hamiltonian is less or equal than E. Finally, h 3 is the phase volume of 
an individual microstate, used to adimensionalize the number of microstates, and must be 
multiplied by iV! in the case of indistinguishable particles. Both the precise value of h and 
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the presence of N are, however, irrelevant to the following discussion. For most systems 
(specially for moderate and large values of N), T is a strongly increasing function of E, 
so r (E — 5E, V, N) <C T (E, V, N) and the number of microstates available to the system 
can be approximated by T (E,V, N), and the probability of each of them by its inverse. 
Substituting this probability in fllap a straightforward ^-generalization of the well-known 
microcanonical entropy equation is obtained: 

S q (E,N,V) = k B \og q T(E,V,N). (4) 

In conventional thermodynamics, temperature is defined as the inverse of the partial 
derivative of entropy with respect to energy: 

T -\~dEj VN -k B V 1 (5) 

where 

These equations, relating entropy and temperature with T and Q, are responsible for the 
differences between the present development and the conventional extensive one. 

The microcanonical average of a general magnitude x (r N ) which depends on the positions 
of the particles is expressed in two equivalent ways: 



(X) = -J wf jx (r N ) e[E-H (r^p*)] df N df (7) 
(X) = j^jx (r N ) S[E-H (t* j?")] dr^df. (8) 

The equivalence is based on the hypothesis that T is a strongly increasing function of E, so 
only the contribution to it from the higher energies, Q5E, must be taken into account, to a 
good approximation. In the remaining of this paper, only the second of these formulae will 
be used. 

In order to obtain a microcanonical formula for any thermodynamic variable, it is neces- 
sary to express it as a function of T or its partial derivatives, like it has already been done 
for temperature. The pressure, isocoric heat capacity and Grunessein parameter are likewise 
straightforwardly written as 
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(9) 
(10) 

(11) 



Other variables often used in tables and experiments, but not so inmediately translated into 
microcanonical language are the isobaric heat capacity, coefficient of thermal expansion and 
isothermal compressibility: 



C„= — 



dE\ 



1 (dV\ 



i 



V 



(12) 



8V 

p,N r \ K ^- L / p,N " \®P/T,N 

This is because of the fact that the constraints more easily applied to a system in the 
laboratory (constant temperature or pressure) are not those of the microcanonical ensemble 
(constant energy or volume). However, the previous variables can be related to 
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using the following identities, obtained from standard relations among derivatives: 



(13) 
(14) 



Kt 



a = kt \ — — cxe-de 



V 



a 



C p = C v [l-—\. 



(15) 

(16) 
(17) 



It is necessary to find a way to express the derivatives of T and Q in such a way that they 
can be calculated directly during a MD simulation. For this purpose, one can observe that 
the integration in ([3]) can be performed the following way due to the form of the hamiltonian 
ineq. ©: 



df dp 



H(r N ,p N ) 



U(r N )<E 



dp 



E E £;<E-U(r") 
i=lj=l 



df 



The inner integral in the right-hand side of this equation is clearly the volume of a 3N- 



3JV 

dimensional sphere of radius R = \j2m\E — U (r N )], which equals 7 3 T\ , R 3N Substituting 
this into the previous equation, eq. ([3]) becomes 



F(E,V,N) 



(2m7r) 2 



9[E-U (f N )] [E-U (f N )] 2 df 



(18) 



Using an almost identical course of reasoning and the formula for the surface (instead of the 
volume) of the n-dimensional sphere, eqs. <Q and ([8]) can be reexpressed as: 



3N 



:'..Y 



<*> = 1 ni / * « 9 I E - 17 I B - 17 ^ **■ (20) 
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These expressions are also valid for the extensive case (since they do not depend on the 
choice of entropy functional as long as the latter preserves the microcanonical distribution) 
and were already provided by Pearson et al. [4] using the Laplace transform method. 

A first practical result can be extracted from these transformations by setting x {r) — 
E — U if) (the kinetic energy expressed as a function of the positions): 



^ = ai^m-i)\ S 0[E - U (r " )] [E - u (fN)] * **■ 

and comparing this to eq. ([TBI , concluding that (E c ) = Combining this with the 

definition of temperature given in (jSJ) an expression for the average kinetic energy as a 
function of the temperature can be formulated: 

(E c ) = —k B TT*-\ (21) 
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which can be considered a g-generalization of the theorem of equipartition of energy. The 
extensive equation is straightforwardly recovered in the q — > 1 limit. Equation (T21T) can be 
used to calculate the microcanonical temperature during a simulation, but it suffers from the 
drawback of having an explicit dependence on T, an integral not usually employed directly 
in the field of MD. T can be calculated, for a given set of parameters, from a numerical 
approximation of integral (Tl8|) . Unfortunately, the complexity of this calculation grows 
exponentially with iV and the integrand is highly discontinuous, which may preclude micro- 
canonical calculations for large systems in the Tsallis formalism. Furthermore, a magnitude 
which depends on the absolute entropy seems unphysical. 

The equation for the pressure can be obtained by substituting ([IS]) in (jH]) (taking into ac- 
count for the derivation that both the limits and the integrand depend on V) and comparing 
with ([6]), which gives: 



By differentiating ([H]) with respect to E and comparing it with the result of substituting 



X = E- 1 = [E-U (f N )] 1 into one arrives at (§§) yiV = fi - l) (E' 1 ), and eq. 



To calculate the next variable of interest, 7, there is only an ingredient lacking according 
to eq. (1111) . the partial derivative of Q with respect to the volume. Using the same method 




(22) 



ffTUj) can be rewritten as 




(23) 



as for the pressure, this parameter is found to be equal to 4r — ^ — 1 
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so 
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No new derivatives are needed to rewrite expression (fTBl for a^: 
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Finally, differentiating (|18|) two times with respect to volume 
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which, after substitution in (fill) , gives the formula for Be'- 
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Partial derivatives with respect to volume can easily be transcribed using gradients as it was 
done in eq. ( |22l ). 

These particular expressions have been chosen so all of them except ( I2ip are formally 
independent of T; hence the only additional difficulty for a nonextensive simulation is the 
calculation of the temperature, which has already been discussed. All the formulae either are 
independent of q, and thus valid also for the extensive case, or recover the SBG expression 
in the g — ^ 1 limit. 

The ideal gas and other simple systems are often used as examples in textbooks because 
their phase-space volumes can be calculated analytically. In particular, for the perfect gas, 
it is possible to take U = 0, so E > and the integrals over real space in the previous 
section become trivial: 
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results which can be substituted in eq. ([5]) to obtain the caloric equation of state: 
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Using (1221) and ( I2TT) . the thermal equation of state can likewise be obtained: 
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It is easy to see that as q — > 1 these equations tend to the well-known extensive ones. 
Important simplifications arise also when N is large. These two limits, however, do not 
commute. More about this can be found in the canonical treatment by Abe in Q|. A differ- 
ent microcanonical treatment has been proposed by Parvan |6(, who includes an extensive 
thermodynamical variable £ = - — 1 to the description of a system in order to recover the 
extensivity in the thermodynamic limit. 



III. FORMULAE OBTAINED USING THE PHYSICAL TEMPERATURE 

The consideration of the difficulty of obtaining T by performing the integration in eq. ffT8j) 
may suggest the idea of setting the system at a given temperature by putting it in contact 
with a thermostat before starting the microcanonical simulation, so both T and (E c ) are 
known and F can be calculated using eq. (12ip . However, this approach fails because the 
condition of thermal equilibrium between two systems does not imply that they are at the 
same temperature defined by eq. (J5j), as discussed by Abe et al. in Ref. [3]. Thus this 
definition can be rejected as unphysical. The correct physical temperature can easily be 
shown to be 

1 + ^^-S q j T = T^T, (30) 



which has as an inmediate consequence, taking equations (jSJ) and ( 12T1) into account, that this 
physical temperature, as a function of phase space integrals, is independent of q, as noted 
by Toral [8|. Thus, this temperature can be seen as the usual derivative of the Boltzmann 
entropy. This has led Gross [9J] to conclude that the Boltzmann entropy is enough to describe 
every hamiltonian system at equilibrium, extensive or not. This conclusion has been disputed 



by Wang 10J, and is not applicable to the aforementioned Parvan formalism 6|. However, 
a detailed discussion is clearly outside the objectives of this paper. 

With this definition, the equipartition theorem recovers its familiar form: 
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In particular, eq. (I3T]) is identical to the one obtained using the SBG entropy. Thus, 
the expressions given by Pearson for the pressure, isocoric heat capacity and Griinessein 
parameter, as well as the auxiliary variables a B and Be, or equivalently the q — > 1 (T — > 
Tphys) limits of the ones given in the previous section, are still valid: 
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The only obstacle to simply adopt all the extensive formulae is the definition of pressure 
as a derivative of entropy. To keep its physical meaning, the new definition of temperature 
requires that pressure is expressed as 
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which means that derivatives involving both T p h ys and p must be changed in order to reflect 
this fact. This affects thermodynamic equations ( JTBT) . (JTBI) and ( fPTj) . Since T p h ys is a function 
of S q and T, it is easy to use the chain rule to relate the physical variables to the unphysical 
ones for which equations have already been developed in the previous section: 
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These magnitudes are not independent of q and ultimately depend on T through Cy, which 
means that knowledge of T is still necessary to completely determine the thermodynamics 
of the system. 



IV. CONCLUSIONS. 



The thermodynamic behaviour of a system described by the Tsallis entropy has been 
studied in the microcanonical ensembles. Two possible definitions of temperature have 
been considered, resulting in two different sets of formulae. These microcanonical formulae 
can be used to obtain thermodynamic information (heat capacities, pressure, dilatation 
coefficients, etc.) from mechanical averages available during a MD simulation. However, if 
the first definition of temperature is used, which amounts to a direct generalization of the 
conventional one, these thermodynamic magnitudes depend on phase-space integrals which 
are not easily calculated. The second definition of temperature is designed to satisfy the 
zeroth law of thermodynamics and addresses this problem but only partially, returning some 
formulae to their extensive (q = 1) form. This is explained by the fact that the physical 
temperature is, for any value of q, a derivative of the SBG entropy, a relation clearly seen 
in the microcanonical formalism. 
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